How Should We Modify the High Energy Interaction Models ? 
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An analysis has been made of the present situation with respect to the high energy hadron-nucleus and nucleus- 
nucleus interaction models as applied to cosmic rays. As is already known, there are inconsistencies in the inter- 
pretation of experimental data on the primary mass composition, which appear when different EAS components 
are used for the analyses, even for the same experiment. In the absence of obvious experimental defects, there is 
a clear need for an improvement to the existing models; we argue that the most promising way is to enlist two 
effects which should be present in nucleus-nucleus collisions but have not been allowed for before. These are: a 
few percent energy transfer into the EAS electromagnetic component due to electron-positron pair production or 
electromagnetic radiation of the quark-gluon plasma and a small slow-down of the cascading process in its initial 
stages associated with the extended lifetime of excited nuclear fragments. The latter process displaces the shower 
maximum deeper into the atmosphere. 



1. Introduction 

One of the possible ways to assess the quality 
of the interaction models is to use them for the 
determination of the primary mass composition 
from the analysis of different cosmic ray compo- 
nents. The ideal case is to get a consistent mass 
composition using the same model and different 
components. As the variety and precision of ex- 
perimental data and the quality of the theoretical 
interaction models are improved, more evidence 
appears that their quality is still not good enough. 
The mass composition derived with the use of 
electromagnetic and muon components is system- 
atically lighter than when hadrons and muons are 
used. Mass composition in the knee derived us- 
ing measurements of the depth of maximum is 
lighter than that from all ground-based measure- 
ments [0. All these facts indicate that the exist- 
ing model needs further improvement. 

2. The evidence from the ground-based 
measurements 

The variety of observed EAS components and 
precise measurements of their characteristics by 
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KASCADE experiment allowed this collaboration 
to derive the primary mass composition using the 
multivariate analysis of data. It was found that 
not only do different methods applied to the same 
set of observables give different results, but that 
there is a systematic difference between the re- 
sults obtained using the same method, but dif- 
ferent components. In Figure 1, taken from 
the mean logarithmic mass (InA) derived from 
different sets of observables is shown vs. the so 
called truncated muon number Njj~ , which is an 
observable adopted as a measure of the primary 
energy, independent of the primary mass. The 
fundamcntsl problem is that all the (InA) values 
should be the same at the same value of NfT , and 
they are not. 

A number of features are of interest and, pre- 
sumably, of some importance: 

(i) the use of the electron size N e results in an 
average lighter composition ( Fig. la ). On the 
other hand, omitting electrons and using just 
muons and hadrons results in a heavier compo- 
sition ( Fig. lb ); 

(ii) considering the mean values from the two sets 
of data ( filled-in circles ), the ratio of ((InA)) 
for(b) to (a) increases smoothly with lgN^ r , i.e. 
the discrepancy between (b) and (a) rises with 
increasing primary energy ( Fig.lc ). This dif- 
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§gure 1 . Mean logarithmic mass (In A) from the analysis 
of different sets of observables vs. IgNj^". The two 
dashed lines indicate alternative predictions of the Single 
Source Model H , with the upper line being our latest and 
'best estimate'. The sets displayed in (b) do not include 
the observable N e . It is seen that omitting N e results in 
a heavier mass composition, (c) The ratio of mean values 
of (In A) from sets (b) and (a). 



ference points to inadequacies in the models used 
for the analysis of the experimental data and pro- 
vides an impetus to correct them. However the 
corrections should not be radical, because the dif- 
ference between the (InA) is not large, typically 
8(lnA) ss 0.2 - 0.4. The mentioned systematic 
difference can be an indication that the energy 
distribution between the different shower compo- 
nents is slightly different from that in the models: 
specifically, the actual mean number of electrons 
iV e in EAS appears to be higher, that of muons 
Nf,. slightly lower and that of hadrons Nh lower 
still than in the models. 

It is well known from the present models that 
for the same primary energy the number of 
muons in nuclei-induced showers is higher than in 



proton- induced ones. On the other hand the num- 
ber of electrons and hadrons in nucleus-induced 
showers observed in the lower half of the atmo- 
sphere is lower than in proton ones. If one has 
an opportunity to measure the primary energy, 
by muons, Cherenkov light or another technique, 
and finds that the shower has a low or a high 
N e ( actually the ratio jf- is important ) the con- 
clusion will be that this shower is initiated by a 
proton or light nucleus. On the contrary, if one 
finds that the shower has a low Nh the conclusion 
will be the opposite, i.e. that the shower is initi- 
ated by a heavy nucleus. This is exactly what is 
observed in the showers at sea level. 

The same conclusion can be drawn from an 
analysis of KASCADE event rates [Q. Both the 
muon and hadron trigger rates, observed by KAS- 
CADE, are lower than expected from the model 
calculations. This discrepancy indicates again 
that the actual numbers of muons and hadrons 
in EAS are lower than in the models, although 
the energy region responsible for the trigger rates 
is lower than that analysed for the mass com- 
position around the knee. The analysis of this 
discrepancy indicates that the needed reduction 
of the number of muons in the model should be 
about 6%, for hadrons it is bigger - about 29%. 
Because muons and hadrons are the products of 
hadronic cascades, it is evident from the energy 
balance that to reduce the energy contained in the 
hadronic cascade one has to increase the energy 
transferred into the electromagnetic cascade. 

3. Numerical estimates 

3.1. Balance between the EAS components 

To evaluate the effect of the proposed change 
of the balance between different cascade compo- 
nents we applied the semi-quantitative analytical 
approach, details of which can be found in j^j. 
The results of the calculation for our basic set 
of parameters are shown in Figure 2a by a full 
line. According to our suggestion we increased 
the energy fraction transferred by nucleons into 
the electromagnetic component K N from 0.20 to 
0.26. The result is shown in Figure 2a by the 
dashed line: the muon energy at sea level de- 
creased by ~ 6%, the hadron energy decreased 
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by ~ 23%, the energy transferred into the elec- 
tromagnetic component increased by ~ 2%. 



LONGITUDINAL DEVELOPMENT OF 1 PEV CASCADES IN THE ATMOSPHERE 
Development of the nucleon (N). pion (tt). muon and electromagnetic components Ce?) 




Figure 2. The longitudinal development of 1 PeV cas- 
cades in the atmosphere: (a) Fractions of the total en- 
ergy carried by nucleons (N), pions (it), muons(/i) and 
transferred into the electromagnetic (e*y) component; (b) 
Electron size of the shower N e . Basic parameters: full 
line: K*f = 0.20, ER = 65 gem" 2 ; dashed line: K!* = 
0.26, ER = 65 gcm~ 2 ; dotted line: = 0.26, ER = 71 
gem -2 . 



3.2. The triangle diagrams 

The balance of the energy contained in the ma- 
jor EAS components is convenient to analyse us- 
ing the so-called 'triangle diagrams' ||. If the 
height of an equilateral triangle is equal to 1, then 
for each point inside this triangle the sum of the 
distances to its sides is equal to 1. If we know the 
energy fractions carried by the electromagnetic 
(<5 e7 ), muon (5^) and hadron (5h) components at 
the observation level, so that <5 e7 + 8^ + Sh = 1, 
then each shower can be presented by a sin- 



gle point inside the triangle. Our basic shower 
( Fig. 2, full line ) is shown in Fugure 3 by a full 
circle. The desired direction for the shift of the 
energy balance in the modified model is shown by 
the straight line arrow in Figure 3b. 



TRIANGLE DIAGRAM FOR THE ENERGY BALANCE OF 1 PEV EAS AT SEA LEVEL 
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Figure 3. Triangle diagram for a 1 PeV shower at sea 
level, (a) Large scale diagram with the inset magnifying 
the indicated region, (b) Small scale part of the full di- 
agram. The straight arrow to the left of the full circle 
indicates the desired modification of the EAS energy bal- 
ance. 



However, despite the increase of the energy 
transferred into the electromagnetic component, 
the 'preserved electromagnetic energy and the 
electron size of the shower at sea level decreased 
by 14% due to its faster development and then 
faster attenuation of the cascade ( Figure 2b ). 
The point in the triangle diagram moved in the 
different direction ( open circle ). The similar re- 
sult is obtained if we increase the total inelasticity 
Kf ot with = ^K^ ot . Therefore the mere in- 
crease of the mean value of cannot give the 
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required result. 

We argue that another effect that should be 
present will bring about the desired effect: the 
slowing down of the development of the cascade 
in its initial stages. For illustration purposes we 
slow down the development of the hadronic and 
electromagnetic cascade by increasing the elon- 
gation rate ER from 65 gcm~ 2 to 71 gem -2 , pre- 
serving K*? — 0.26. The result is shown in Figure 
2 by dotted lines. The direction and the mag- 
nitude of the changes is now correct. Thus we 
conclude that the increase of the energy trans- 
ferred into the electromagnetic component com- 
bined with the slowing down of the development of 
cascades in their initial stages is the most realis- 
tic way to improve the particle interaction model 
and to achieve a consistent estimate of the pri- 
mary mass composition. 

4. Theoretical arguments 

Besides all these arguments, which are purely 
phcnomcnological, there are also theoretical ar- 
guments which lend support to the phenomeno- 
logical consideration. Nearly all of them are 
related to processes which appear in nucleus- 
nucleus ( AA ) interactions. Their details are 
given in [g| and here we just enumerate them. 

In AA collisions with a small impact parame- 
ters ( central collisions ) one has to expect the 
production of e + e~-pairs including the multiple 
pair production. This process provides an ad- 
ditional energy transfer into the electromagnetic 
component and it was not taken into account in 
the present models. An additional energy trans- 
fer into the electromagnetic component can arise 
also from an excess of direct photons, which has 
been predicted theoretically as a signature of the 
quark-gluon plasma and is now observed in AA- 
collisions both at low and high transverse momen- 
tum. 

In AA collisions with a large impact parame- 
ters ( peripheral collisions ) a projectile nucleus 
fragments into few pieces of different mass. Some 
of them are excited and after de-excitation give 
rise to MeV gamma-ray lines, observed from 'dis- 
crete' sources and the interstellar medium. The 
lifetime of the excited fragments varies from a 'nu- 



clear' time ~ 10~ 23 sec to millions of years. For 
AA-interactions at PeV energies both the lifetime 
before the de-excitation and the energy of emit- 
ted gamma-quanta are extended by the factor of 
10 5 - 10 6 due to relativistic effects. As a con- 
sequence one can expect an additional sub-PeV 
electromagnetic cascade to be initiated a few hun- 
dred meters below the point of the first interac- 
tion. This effect will slow down the development 
of the electromagnetic cascade and shift its max- 
imum. 

5. Conclusions 

The inconsistencies in the interpretation of the 
experimental data on the primary mass compo- 
sition, obtained when different EAS components 
are used for the analysis, indicate the need for 
some improvement to the models which were used 
hitherto. We propose that the most promising- 
way is to introduce an additional ( a few per- 
cent ) energy transfer into the EAS electromag- 
netic component combined with a slowing down 
of the cascade development in its initial stages, 
which is followed by a small ( 20-30 gcm~ 2 ) shift 
of the shower maximum into the deeper atmo- 
sphere. The most likely processes which can be 
responsible for such changes are those which oc- 
cur in AA-collisions and they should indeed be 
present at some level. The importance of these 
processes is expected to grow with energy and of- 
fers the hope of resolving some controversies at 
very high energies. 
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